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Abstract: The holographic description of Pomeron exchange in a strongly-coupled gauge 
theory with an AdS dual is extended to the case of two to three scattering. We study 
the production event of a central particle via hadron-hadron scattering in the double Regge 
kinematic regime of large center-of-momentum energy and fixed momentum transfer. The 
computation reduces to the overlap of a holographic wave function for the central particle with 
a source function for the Pomerons. The formalism is applied to scalar glueball production 
and the resulting amplitude is studied in various kinematic limits. 
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1. Introduction 

Regge phenomena have been studied, both in experiments and in theoretical contexts, for 
several decades. Experiments investigating hadron scattering, and associated theoretical at- 
tempts to understand the data, led in the 1960s to a number of important developments that 
continue to play a role in current research. The general Regge approach to writing amplitudes 
as complex- valued functions of Mandelstam variables, and investigating them using Mellin 
transforms, had some general success in characterizing observed amplitudes that appeared 
at the time to be inconsistent with quantum field theory. The development of "dual reso- 
nance" models that captured various features of scattering amplitudes, including s-t channel 
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duality, towers of resonances lying on Regge trajectories a(m^) for m? > 0, and near-forward 
scattering amplitudes proportional to s"^*^ for t <0, were the first steps along the road to a 
consistent string theory. Of course, hadronic scattering experiments and standard string the- 
ory diverged in large-angle scattering (where the amplitudes for the former fall as powers of 
s, while those of the latter fall exponentially) and in deep inelastic scattering (where Bjorken 
scaling suggested hard, weakly-interacting substructure). The theory of QCD supplanted 
string theory as the leading model for strong interactions, and its success is spectacular. Fur- 
ther, it was soon learned that consistent string theories needed to be defined with massless 
spin-two particles, which QCD of course lacks. Yet the fact that Regge phenomena are evident 
in the data, in regimes which are not amenable to study using perturbative QCD, and the 
successes of the early phenomenogical (and inconsistent) string theories, have left lingering 
questions about the relationship between string theory and the physics of hadrons. 

A better understanding of the relationship between QCD and string theory has now 
emerged from the discovery of an apparent (but still unproven) duality between non-abelian 
gauge theory and string theory in curved spaces. This duality conjecture [1] contends that 
four-dimensional gauge theories which are asymptotically conformal are identical to ten- 
dimensional string theories propagating on spaces that are asymptotically five-dimensional 
anti-de Sitter space (AdS^) times a compact five-manifold W. The conjecture has provided 
a setting in which the successes and failures of the old string theories could be substantially 
clarified. The inability of four-dimensional fiat-space string theories to match general features 
of high-energy QCD amplitudes was shown to be rectified in the scattering of strings propa- 
gating on the appropriate curved spaces; exponential fall-off with s in large-angle scattering 
is replaced with power laws [2]. Bjorken scaling does not hold, but its failure is replaced [3] 
with a more generalized scaling discussed by Kogut and Susskind in the context of conformal 
field theory [4]. The Regge phenomena in string theory amplitudes manage to reproduce 
diverse phenomena of gauge theory in various Regge- like regimes [2,5]. In high-energy scat- 
tering at small momentum transfer \t\, where experiments show the classic Regge forward 
peak, the string theory computation reduces approximately to four-dimensional flat-space 
string theory, and the Regge peak in the string theory is transferred directly to the gauge 
theory. At larger negative t, however, and in some other regimes as well, the fifth dimension 
of the asymptotically- AdSs space begins to play a role, washing out the stringy Regge behav- 
ior, and instead reproducing phenomena that mirror the amplitudes calculated by Balitsky, 
Fadin, Kuraev and Lipatov (BFKL) [6] in gauge theory, as well as other phenomena. At the 
same time, the string theory also provides a discrete set of states lying on a set of Regge tra- 
jectories. Thus many of the essential features of gauge theory, especially those which are not 
found in purely perturbative QCD, have been shown to be reproduced in this "gauge/string 
correspondence." 

In this paper, we continue the process of learning about Regge physics in contexts that 
are challenging to understand fully in QCD. All previous studies of Regge physics in the 
gauge/string correspondence have been of 2 ^ 2 scattering, and here we will attempt to 
extend this to important 2^3 phenomena in the so-called "double diffractive" or "double 
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Reggc" regime. In this regime, one considers AB — >■ ABX, where A, B are hadrons that 
are scattered at small angles without internal excitation or other disruption, and X is an 
object produced with a rapidity that lies well between the rapidities of the outgoing A and 
B. Interesting examples in this class that are more or less addressable experimentally include 
cases where A is a proton, B \s a proton or antiproton, and X is a glucball, a quarkonium 
state, or a Higgs boson. Diffractive quarkonium production has been observed at Fermilab [7]. 
Higgs production in double diffractive pp scattering is an important controversial topic, as 
experiments are being planned to search for it at the LHC amid heated discussions about 
the expected cross-section (see [8] for a brief review). We will however focus on something 
akin to pp ^ pp + glueball, as this is technically the easiest problem. Our methods will 
require some minor generalization for the study of quarkonium production, and at least one 
additional technical advance for Higgs production. 

We should emphasize that our goals are strictly limited, yet ambitious. On the one hand, 
none of these processes in real- world QCD can be reliably computed using the gauge/string 
correspondence. The correspondence allows detailed study of a certain class of gauge theories, 
but unfortunately QCD is not among them. No well-controlled study of QCD dynamics, al- 
lowing for a computation of cross-sections for comparison with experimental data, is currently 
possible. While some surprising success has been achieved modeling the spectrum of QCD, 
as well as some hadron couplings, using five-dimensional effective field theories [9, 10], it is 
quite another matter to compute scattering amplitudes in which stringy phenomena play an 
essential role. We know that some aspects of QCD in scattering are not true in gauge theories 
for which stringy computations are reliable. For instance, Bjorken scaling is badly violated, 
and Regge exponents differ by something of order unity from those at weak coupling. 

Our approach here is to treat those gauge theories for which the gauge/string correspon- 
dence can be easily applied as toy models for QCD. We seek to identify phenomena which 
are universal or quasi-universal in these toy models, and in turn, to understand the degree 
to which they may apply also in all confining gauge theories, including QCD. Though the 
numerical details of our calculations will not match QCD, we expect our toy model and 
QCD to share key qualitative and semi-quantitative features. Moreover, the similarities and 
differences between our model and QCD should be physically comprehensible. 

The simplest possible use of our toy models would be a direct one: it is possible that 
the behavior of our amplitudes as a function of the kinematic variables may be similar to 
that of QCD. A more subtle but potentially more important use of these toy models is at the 
level of general methodology and conceptual understanding. Our goals here are limited to 
the latter, and our results should be viewed as exploratory, especially as data on the process 
computed below is limited. We hope this paper will be a useful step in the direction of 
allowing computation of diffractive processes for which data is more easily accessible. 

2. Preliminciries 

A gauge/string duality is a conjectured equivalence between a string theory in an asymptoti- 
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cally hyperbolic geometry and a gauge theory in a fewer number of dimensions. The original 
and best established example of this duality maps type IIB string theory in an AdS^ x 
background, where AdS^ is five-dimensional anti-de Sitter space and is a five-dimensional 
sphere, to the maximally supersymmetric SU{N) Yang-Mills theory in 3-1-1 dimensions. This 
supersymmctric field theory is scale invariant and is thus an example of a conformal field the- 
ory (CFT). As the hard- wall model we consider in this paper is only a slight generalization of 
this original version of the duality, it is worthwhile to review the details of the AdS^ /CFT4 
correspondence. 

The original gauge/string duality can be motivated by considering the effects of placing 
a stack of A'^ coincident D3-branes in flat ten dimensional space. The D3-branes are 3-1-1 
dimensional surfaces on which the open strings of type IIB string theory may end. On the 
one hand, a low-energy description of a stack of N D3-branes is = 4 supersymmetric 
Yang-Mills theory, loosely motivated by the fact that the open strings which interconnect 
the D3-branes can be reinterpreted as gluons. On the other hand, the D3-branes are massive 
objects and in the large N limit substantially warp spacetime. "Close" to the D3-branes, the 
geometry approaches AdS^ x S^. 

Given this correspondence, for which there is now an enormous amount of evidence but 
no proof, there exists a dictionary mapping stringy quantities to gauge theory quantities. 
The Yang-Mills coupling constant gyj^^ maps to 'ings where gs is the string coupling constant 
which measures the probability for strings to split. The 't Hooft coupling A = gyM^ maps 
to /a'"^ where R is the radius of curvature of both AdS^ and and 1/2^ a' is the string 
tension. One usefulness of the duality lies in the scaling limit N 00 while keeping A large 
and fixed. By taking A^ 00, the string splitting amplitude is suppressed and strings may be 
treated classically. Taking A large corresponds to keeping the curvature scale of the geometry 
very large compared to the string scale, in which limit low-energy processes in string theory 
are well approximated by supergravity. Thus strongly interacting physics on the gauge theory 
side gets mapped to classical general relativity. 

In this paper, we will be interested in very high energy scattering processes for which su- 
pergravity is not enough, although it provides a useful point of departure. While the original 
gauge/string dualities have only CFT duals, certain generalizations exhibit renormalization 
group flow and confinement. The hard-wall model we consider below is a trivial generaliza- 
tion of the ^dS'5/CFT4 correspondence where a confinement scale is introduced by hand by 
removing a portion of the interior of AdS^. We will consider high energy 2—^3 scattering 
processes in this cut-off geometry. 

Regge phenomena in gauge theories well-described by the gauge/string correspondence 
were studied by a number of authors. The methods used in this paper appeared first in the 
small-x region of deep inelastic scattering [3] and were fleshed out much more fully in [5]. 
There, following ideas of ref. [3], the Regge limit of 2 — >^ 2 scattering was studied. In a regime 
where N is taken very large and A and s are taken very large in a correlated way, with t flxed, 
the scattering amplitudes are dominated by single-Pomeron exchange. The Pomeron is a 
coherent color-singlet object built from gluons whose properties are universal; it is the object 
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which is exchanged by any pair of hadrons that scatter at high energy and large impact 
parameter. In the dual string theory, the Pomeron is not the graviton but the graviton's 
Regge trajectory. 

While string theory in flat space in the Regge limit exhibits Regge scaling s°'^^\ in the 
gauge/string duality context this flat-space string amplitude must be corrected by the warping 
of the metric and flnally convolved with supergravity wave functions representing the scattered 
hadrons. Moreover, as was argued in ref. [5], when ln(s/A^) (A of order the confinement scale) 
becomes large compared to \/A, the local form of the string amplitude becomes inappropriate, 
and the local Mandelstam variable t must be promoted to a differential operator that acts on 
these supergravity wave functions. This in turn leads to diffusive behavior of the scattering 
strings. This diffusion is the strong-coupling analogue of the diffusive behavior seen at small 
A in the calculations of BFKL [6]. In gauge/string duality, the radial direction of AdS^ has 
a dual interpretation as an energy scale in the field theory, and the diffusion happens in this 
radial direction as well as in the directions tranverse to the scattering direction. In BFKL, 
the role of the radial dimension is played by a transverse momentum variable circulating in 
the ladder graphs of the Regge resummation. 

In this paper, we apply the methods of ref. [5] to a more elaborate problem. We study 
2^3 scattering, in a regime where two Pomerons are emitted by the initial-state particles 
and fuse to make a third particle. Here we will consider this particle to be a glueball: 
a normalizable state in the AdS^ space. Experimentally one might also be interested in 
quarkonium states, or a Higgs boson, but we will not consider these in this paper. 

We begin in section ^ by reviewing the kinematics for these 2 — > 3 processes. Instead of 
two independent Mandelstam variables s and t, five-point amplitudes involve the variables s, 
si, S2, ti, and t2- The double diffractive limit, or double Regge limit, we consider consists of 
taking s and Si very large compared to while keeping the ti small. In double diffractive 
scattering in the center-of-momentum frame, the hadrons scatter by very small angles. Since 
we will explore this 2 — > 3 process in the hard-wall model, we review details of this model in 
section ^. We also present the glueball and hadron wave functions that will be used in the 
calculations. 

The next step is to construct the five point, fiat space, string theory amplitude which we 
do in appendix While the four point amplitude, or Virasoro-Shapiro amplitude, is well 
known and relatively simple, the corresponding closed string five point amplitude is vastly 
more complicated and less well known. In general, using the techniques of [11], the five 
point amplitude can be expressed as a quadratic polynomial in generalized hyper geometric 
functions 3-F2- Fortunately, we only need this amplitude in the double diffractive limit, and 
our discussion is taken in large part from [12]. 

In section ^, given the fiat space amplitude in the double diffractive limit, we use the 
procedure outlined in [5] to convert this fiat space amplitude into a curved space amplitude. 
The amplitude may be expressed as an integral of the glueball wave function over a source 
function TZ which is a property of the two fusing Pomerons. Many properties of this 2 — > 3 
process can be evaluated independently of the produced fifth particle by studying TZ. 
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In section g we pause to explain the constraints placed on the various parameters of the 
scattering process by the kinematics. 

We evaluate this double difFractive scattering amplitude in various regimes of tj and Si in 
section |^. We begin by considering ti = for large values of ln(s/A^) and ln(si/A^). In this 
case, we find that the scattering amplitude, an even function of rapidity, y, is a function only 
of the combination y/^/\. In other words, at large A, the amplitude is nearly independent of 
y. This is consistent with the corresponding dual gravity amplitude. 

Next, we consider arbitrary tj and ln(s/A^), ln(sj/A^) ^ 1. For large ti we find power- 
law behavior fall-off with tj, without any sign of the exponential Regge peak in the forward 
region. The Regge peak is only observed for smaller values of U and ln(sj/A^), as a transient 
rather than an asymptotic phenomenon. This is roughly consistent with the hard scattering 
results of [2]. 

Finally, we present a concluding discussion in section ^ 



3. Double diffractive kinematics of 2 ^ 3 

In this paper, we calculate 2^3 scattering amplitudes in the double diffractive limit. We 
have in mind a process similar to pp ^ pp + glueball. The incoming hadrons are defiected by 
very small angles from their original trajectories and produce a glueball via double Pomeron 
exchange. 



->- 
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Figure 1: Momentum flow diagram for double diffractive scattering. The hadrons have momenta 
ki, . . . ,k/^ and the glueball has momentum k^, . Pomerons are represented by the internal lines in the 
ti and t2 channels. 

Working in — h -|— |- signature, we will use generalized Mandelstam variables 



s 

si 

h 

t2 



-{ki-k2)\ 
-{ks-k^f, 



(3.1) 
(3.2) 
(3.3) 
(3.4) 
(3.5) 
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and the five masses mf = —kf. Our conventions are shown in Figure |T|. We will assume that 
the energies are large enough to neglect the external hadron masses. The glueball has mass 
m5 = m. 

The double diffr active limit we consider is the double Regge limit [13] defined as 

s,si,S2— >oo with s^si,S2 and ti,t2, fixed. (3-6) 

s 

In this limit, and going to the center-of-momentum frame, it is helpful to parametrize the 
momenta in the following way: 

^2 = (^,0,0,^), k3 = {E,0,0,-E), 
ki = + xks + h±, h = + + h±, 
^5 = (1 - e - X)k2 + (1 - X - Oks - {kiL + A;4±). (3.7) 

The momenta and A;4_l are defined to be orthogonal to k2 and ^3. The eight scalar 
parameters introduced above must be expressable in terms of Lorentz invariants and m. 
Putting momenta on-shell and solving algebraically gives 

E = ^/2, x = -h/s, x = -t2/s, 

C = 1 + (ti - S2)/s, e = l + (i2-si)/s, 
klj_ = -h - ti{ti - S2)/s, klj_ = -t2 - t2{t2 - si)/s, 
2fci± • fc4± = ^ - m2 + 2*1*2 + t2 (1 - f ) + ti (1 - f ) . (3.8) 

The conditions on the Mandelstam variables in the double diffractive limit translate into the 
relations ^, ^ ~ 1 and X)X ^ 1- These relations in turn can be understood physically as the 
fact that the hadrons are deflected by very small angles in this limit. 
It is also convenient to define 

mi =m2 + (fci_L +A;4_l)^ (3.9) 

This is the effective mass of the glueball, as viewed in the x~^-x~ light-cone plane. In terms 
of Mandelstam variables, 

2 •5l'S2 f ^ h — t2\ f ^ , h — t2 \ _ S1S2 



mi = ^[l- ^ ^ 1 + ^ ^ « ^ (3.10) 

in the double diffractive limit. 

The rapidity y of the glueball is defined by tanh y = k^z /E^. In terms of the Mandelstam 
variables, the rapidity can be expressed as 

2 \E5-k5J 2 \si+ti-t2j 2 \sj 

In the center-of-momentum frame, for a massless particle, the rapidity is related to the polar 
angle (relative to the beam axis) via = cot(0/2): large \y\ corresponds to ~ or tt. 
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4. Hard-wall model 



The hadronic wave functions that are needed in our computations depend on how confinement 
is incorporated into the gauge/string framework. In this paper, we choose to use the hard-wall 
model. ^ The model imposes a sharp cutoff on the AdS^ radial coordinate z at some scale zq, 



E? 

ds^ = —{rjfj^iydx^dx" + dz^) + ds^, Q < z < zq. (4.1) 



The line element dsy^ gives the metric on the five-dimensional transverse space, which in 
simple examples is just . The cutoff establishes a mass gap A = I/zq- To a four-dimensional 
observer, z ^ is the UV, while z ^ zq is the IR. 

Scalar hadron wave functions are solutions to the ten-dimensional Klein-Gordon equation 
using the metric (4.1). By separation of variables the solutions can be written as e^^'^4'{z)f{6), 



where k is the gauge theory four-momentum. For simplicity, we take / to be a constant mode 
on the compact space W. The function / plays no role in our calculations. When we introduce 
diffusion kernels, it will be convenient to use the dimensionless coordinate u = ln(zo/z), 
< n < oo. To a four-dimensional observer, u — > is the IR and u — > oo is the UV. 

The radial- "Kaluza-Klein" normalizable mode corresponding to a hadron is given by [15] 

^^"^^ ~ AVm^R'/' JA-2(m/A) ' ^^-'^ 

where Voliy is the volume of the compact space W, A is the conformal dimension of the dual 
gauge theory operator, and m is the four-dimensional mass given by m/A = CA-3;n- Ck;n 
denotes the nth zero of the Bessel function Jk{x)? 

In Figure |^, the hadron with momentum represents the glueball. In non-Abelian 
gauge theory, a scalar glueball can be created by the operator O = TvF^uF^^^ . According to 
the AdS'/CFT correspondence, this operator is dual to a massless closed string dilaton state 
propagating in the ten-dimensional bulk whose AdS radial profile is given by the normalizable 
mode (U) with A = 4. 

The external hadrons with momenta fei , . . . , /C4 in Figure |l] are all scalar normalizable 
modes (f)Q of the form ( [4.2[) with A = Aq and mass mo- Although they are not baryons, and 
are scalars rather than fermions, they are reasonable surrogates for protons, as their profile 
in the radial AdS direction is appropriate for an object built out of a small number of valence 
partons (and created by an operator of small twist [2]). 



^The hard-wall model is not a fully consistent SUGRA theory as it does not satisfy the supergravity 
equations of motion. However, past calculations have shown that the model is useful for capturing model- 
independent behavior of confining gauge theories [2,3,5,9,10,14]. 

^We use the generalized Neumann boundary conditions du (^e~^^~'^^'^(l){u)^ = at « = employed in [15], 
which imply (f)' {Q) = (4 — A)(f){Q). These reduce to ordinary Neumann boundary conditions only in the case 
A = 4. The normalization is fixed by the condition VolM' i?^A^ Jp°° du e^"(/)(it)^ — 1. See the appendix of ref. [2] 
for a derivation of this normalization condition. 
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5. S-matrix 

We want to compute the S-matrix, S, for hadron scattering in four dimensions. In string 
theory, at leading order in 1/N and 1/\/A, this is given roughly by a path integral over a 
spherical worldsheet embedded in the cut-off AdS^ x W space, with appropriate vertex oper- 
ators representing the external hadron states. For scattering processes in d flat dimensions, 
we also define the amplitude 7^ by 

5 = i(27r)'^5'^(^A:,)Trf. (5.1) 

Using a prescription given in [5] we determine S by integrating Tiq (the scattering amplitude 
for closed strings in ten-dimensional Minkowski space) with the product of hadron wave 
functions over all coordinates in the cut-off AdS^ x W, 

S= [ d^xdz [ d^9V^Tioiki,...,h)f{e'^^-''(l)iiz)- (5.2) 

JAdSs Jw 

The factor \/—G is the square root of the determinant of the metric on AdS^ x W. It is 
important to distinguish the two sets of momenta that arise in (|5.2| ). The ki are the four- 
momenta of scalar hadrons in the dual gauge theory defined on the boundary of AdS^, while 
the ki can be interpreted as (center-of-mass) ten-momenta of closed strings propagating in 
the AdS^ X W background at a given value of z. The relationship between ki and ki will be 
discussed later in this section. 

Our first task is to compute the tree-level string amplitude Tiq for five external states. 
Rather than use the closed superstring, we calculate in appendix |A| an equivalent amplitude 

for the closed bosonic string. In the double Regge limit, where exchange of the leading 
Regge trajectory in the ti and t2 channels dominates, we expect bosonic and superstring 
amplitudes to be essentially identical, up to overall constants and subleading effects which are 
not important for us here. Moreover, in this limit the amplitude will not depend sensitively on 
the precise vertex operators used to represent our hadrons on the external legs. In particular, 
only their energy-momentum tensor is important, as a source for the emitted Pomeron. Thus 
there is no cost in using tachyons as the initial state vertex operators in the 2 — > 3 scattering 
process. We must be more precise about the fifth particle (the "glueball") whose coupling 
to gravitons, and therefore to two Pomerons, contributes nontrivial kinematic factors to the 
amplitude. We use the dilaton as the particle produced by the colliding Pomerons; this 
corresponds to producing a scalar glueball in the gauge theory. 

Thus, we calculate the bosonic closed string amplitude ^or four tachyons and one 
dilaton in the double Regge limit, with the dilaton singled out as the produced particle. 
In appendix it is shown that 726 can be written as a power series in the dimensionless 
parameter a'm'j_. Here fh'^, the ten-dimensional generalization of m^, is approximately a ratio 
of ten-dimensional Mandelstam invariants, S1S2/S, in the double Regge limit. As we will see 
in section |6|, the physics forces us into a regime where this parameter is much less than 1. 
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The leading term in the series for is given in eq. ( A.26 ). We modify the normahzation 
of ( |A.26| ) so that it may be interpreted as a ten-dimensional scattering amplitude, taking 
gc ~ ct'^gs- To leading order in a'm^ and up to a numerical factor, 



Tio ~ a'^gl 



(5.3) 



'i-K/2 



2 + Q'tl/2 



a s 



+ 



2+a't2/2 



,^m/2'^_^ 



t/2« ^1 
4 



a'(t2-ti)/2 



n(a'ti/4,a'(t2 -ti)/4) 



a'(ti-t2)/2 



n(a't2/4,a'(ti -t2)/4) 



where 



Yl{x,5) = {x + 5)- 



n-5) 



T{2 + x) T{l + 6y ^^-^^ 
It was shown in [5], following [2,3], that a string amplitude in flat spacetime can be 
carried over in some circumstances to an amplitude in a weakly-curved spacetime. To do this 
one must relate the ten-dimensional Lorentz invariants in the string theory: s, s'i,S'2;ii)i2 5 
to the four-dimensional invariants in the gauge theory: s, si, S2, ti, ^2- We define an effective 
string length-squared as a' divided by the AdS warp factor,^ 



a 



» = a'z^/R 



eff 



(5.5) 



Each kinematic quantity s,Si,ti in the string theory should be understood as standing in 
for a differential operator, a Laplacian acting on the appropriate states. The eigenvalue of 
this operator in flat space would be the usual kinematic quantity in the string theory, but 
in curved space the Laplacian is not trivial. For example, the quantity a'ti is a differential 
operator 

This operator, the covariant spin-2 Laplacian appropriate for the Pomeron being exchanged 
in the ti channel [3,5], acts on states 1 and 2; more precisely, it acts on the product of their 
wave functions. Similar expressions apply for t2, ^ and but s, si and S2 are so large 
compared to the remainder of the differential operator that we may approximate 

1 



d d 2 

z—z- h tlZ 



+ 0(1/A) . 



(5.6) 



asi 



2 _ ' / \ 

siz = a^g{z)si . 



(5.7) 



and similarly for S2 and s. In summary, we take 



a s 



a Si 
a'ti 



a 



eff 



Z Si, 



a Vp.. 



(5.8) 



This combination arises naturally in the worldsheet path integral. Separating bosonic fields into constant 
zero modes plus stringy fluctuations, e.g., Z{a^,a^) = z + Z' (a^ , a^), introduces an overall factor of R^/a'z'^ 
in front of the worldsheet action. 
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It is essential to retain the full differential operator in the last expresssion because a'ti appears 
in the exponent of a'J, and s is taken to be exponentially large in \/A. 

We can study the properties of a'Vp. by using coordinates in which the operator becomes 
the Hamiltonian for a quantum mechanics problem, in which an analogue particle moves in 
a potential. Let u = \n{zQ/z) and define 

92 



Hi = -\/Aa'V|. 



where, for the AdS metric with z = zqc the effective potential takes the form 

V^{u) 



4 - z^Ue-^^ 



(5.9) 



(5.10) 



For large negative ti the potential grows exponentially as u decreases. In a confining theory, 
where confinement physics modifies the infrared (small u, large z) region, this exponential 
potential provides an infrared cutoff that screens the details of the physics of confinement. 
As tj — > this screen is removed and the precise nature of confinement comes into play. 

Within the simplistic but useful hard-wall model, the potential takes the form ( ^.1[)| ) 
all the way to u = 0, where the space ends, with an appropriate boundary condition. In 
this case the quantum mechanics problem is completely solved for any ti. A complete set of 
eigenfunctions tp,^{ti,u) with eigenvalues E = 4 + u'^ is given in appendix ^ 

In many circumstances it is useful to write the amplitude explicitly in terms of these 
eigenfunctions. For an arbitrary functional F we use completeness of the eigenfunctions to 
write 



du'(f>i{u')(t)2{u') 



(5.11) 



di^F 



'il)y{ti,u)il}l{ti,u). 



Vp^ is defined to act only on the product (t)i4>2- Likewise, Vp^ acts only on 03(/>4. 

Using (^.11|) it is straightforward to express (|5.2| ) explicitly in the hard-wall model as 

S = {2Tif6^ Ma^gl YolwR^ 

X / du—^(j)^{u) I du (t>i{u')(t>2{u ) / du"(l)3{u")(f>4{u")D{u,u ,u'' 
Jo ^0 -'0 ^0 



(5.12) 



where 



D{u, u , u") 
Dq{u,u',u") 



-W2 «eff('")g ^ 



7-1 /' / /I 
L>q[U, u , u 



(5.13) 



diyTp^{ti,u)ipl{ti,u) I dujip^{t2,u)ip2^{t2,u") 



4 + 1^'^ -u? 



4 + 00"^ -v' 
4^/A ' 4VA 



(5.14) 
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Here is a normalization factor which, among other things, will correct for the fact that we 
have used bosonic strings in place of superstrings to formulate the S-matrix. The variables 



appearing in the exponentials of ( 5.14 ) are 

^ [ln(a;g(z)V4)-m/2], 



2^/X 



n 



2VX 



[lnia',siz)sj4:)-i7T/2] . 



(5.15) 



These are all functions of z, but vary slowly with z. They are analogous to diffusion times [5]. 
Since the integrations in ( p.l2| ) are over regions in which z is not exponentially small in 
\/A, this variation is subleading, as long as we choose to evaluate a'^^^iz) at a sensible place, 
where no large logarithms arise. It is natural to evaluate ( 5.15| ) where the integrand of ( ^.12 ) 
is peaked; we will refer to this value of z as Zscatt ■ 

In the double Regge limit, the s-type 
Mandelstam variables are exponentially large in \/A, so the diffusion times have positive real 
parts. 

The interpretation of S given by ( ^.12| ) is straightforward (reading right to left). We 
convolve pairs of hadron wave functions in the ti and t2 channels with the diffusion kernel, 
and take an overlap of the result with the glueball wave function. Indeed, it will be useful to 
define the quantity 



g4M roo roo 

Tl{u) = / du (j)i{u')(j)2{u') / du" (f>3{u")(t)4{u")D{u,u ,u") . 
^0 -^0 Jo 

In terms of TZ, the S-matrix is simply 

' / du (j)^{u)lZ{u) . 
Jo 



(27r)^5^ (J2 h) Ma'^glYolwR^ 



(5.16) 



(5.17) 



Thus the function TZ, which depends on the specific external hadrons chosen and on the 
kinematics, but not on the produced glueball, may be interpreted as a "source" for the 
glueball state. We will refer to it here as the "double Pomeron source function." 



In ( 5.14 ) we may assume that the size of r is sufficient to guarantee that only the eigen- 
values i/,uj <^ A^/^ are important for the evaluation of the integral. Thus, the arguments of 
n are small and we can expand 11, 



lim n(x, 6) 

x,S^O 



1 

6' 



(5.18) 



Rearranging yields 



Do{u,u',u") ^ -4VA / diyip^{ti,u)^p;{ti,u') / duji^^{t2,u)i!*^{t2,u") 



where 



X e 



-T + Ti + T2 



1/2 — UJ^ 



1/2+4 +4: 



2VX 



[ln(aeff(2;scatt)?TT'i/4) - ivr/2] 



(5.19) 



(5.20) 
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6. Discussion of Parameters and Limits 



Our version of the double diffractive process pp ^ pp + glueball is fully described by the 
choice of six parameters N, A, Aq, mo, A, m, and the five kinematic variables s, si, S2, ti, t2- 

In this work the parameters are constrained as follows. The number of colors A'^ must be 
large to ensure that the S-matrix is dominated by the lowest genus worldsheet. The 't Hooft 
coupling A must be large to be consistent with our calculation of closed strings propagating 
in an AdS background with a large radius of curvature. For the glueball, A = 4 and ni/A 
can be any zero of the Bessel Ji function. 

There are also important constraints on the kinematic variables. In the physical region, 
the momenta transfer-squared ti and t2 are negative semidefinite in any scattering process. 
For Regge behavior to be relevant, we will need |si|, |s2| to be very large compared to A^; 
this requires the parameters r, ri , T2 to be large compared to 1/ \/A. 

One apparent assumption constraining the kinematics arises in the calculation of the 
amplitude 726 1 where we assumed a'm^ <C 1 in order to keep just the first term in a power 
series solution. According to ( |5.8D , in a warped metric this becomes the condition 

aMml = « 1. (6.1) 

This would obviously be satisfied for all z if we were to impose /A^ ^ \/A, or equivalently, 
S1S2/S ^ A^-v/A- However, this condition would be much stronger than necessary. In fact, 
the physical process itself assures that a[,g(z)m^ is never larger than one, without additional 
assumptions. 

The logic is the following. If we produce a glueball with mass of order A, then itself 
can only be large if either —ti ^ A^ or —t2 ^ A^. This follows from the definition (|3.9D , 
which, along with the relations ti ~ ^2 ~ ~^4±^ assures that 

< mi < m^ + 4max(|ti|,|t2|) • (6.2) 

Without loss of generality, let us assume that > |t2|- But if —ti ^ A^, then it serves as 
an infrared cutoff on the dynamics, causing TZ{z) to have support only at regions of z which 
are small compared to l/^/—tl. In our computation this arises from the potential ( ^.101 ), 
which develops a — fie~^" oc —tiz'^ barrier at large z, forcing the physics to smaller z. In 
particular, all the eigenfunctions ipu from which TZ{z) is built have exponentially falling tails 
for z > Therefore, if we take m^ large by making large, we find 

z'^ 1 ?7Z^ 4 

= < < _ « 1 (6,3) 

where we have used eq. ( |6.2[ ). Thus in our calculations we do not need to impose eq. ( |6.1D , 
because it follows automatically from the dynamics. 

One corollary of this discussion is that t_|_ <C 1. Clearly, from the definition of Zscatt 
below ( 5.15|) , it cannot be that Tl{zscatt) is exponentially small, except possibly in regions 
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that make exponentially small contributions to amplitudes. Therefore it follows that the 
above constraint on a'm\ applies. From its definition ( |5.20 ), combined with (|6.3|), t±_ is 
therefore always of order 1 / \/A ^ 1 without further assumption. 



7. Calculation in various kinematic regimes 

Our goal in this section is to compute the double Pomeron source function lZ{u). We do 
this first for the limiting case of vanishing momentum transfers, where we will observe a 
minimal rapidity dependence of the amplitude, and infrared insensitivity for sufficiently large 
s,si,S2. We then consider the case of negative momentum transfers, first observing the 
absence of classic Regge behavior at large Tj, and then identifying it at smaller Tj as a transient 
phenomenon. We also find power-law fall-off at large Tj and large U, analogous to what was 
seen in [2]. 

7.1 Zero ti,t2 and large Tj 

We begin by calculating the double Pomeron source function for ti = ^2 = 0. In this regime 
no infrared cutoff protects the scattering process from the details of confinement, and we do 
not expect the detailed results from the hard-wall model to apply generally to all theories. 
Nevertheless, we may hope for universal behavior in some limited settings. We will find that 
the source function is nearly flat in rapidity, for reasons which should apply to all theories at 
large 't Hooft coupling. We will also see that the source function becomes less sensitive to the 
infrared as the center-of-momentum energy increases. While this would be very interesting if 
it were universal, we unfortunately see no reason why this should be the case. 

In order to find IZ we must first calculate the diffusion kernel according to eq. ( ^.191 ). 
This formula calls for the eigenfunctions of the hard-wall model Hamiltonian at ti = t2 = 0, 
which are given explicitly in the appendix by eq. ( |B.7| ). Upon simplifying, one finds that the 
kernel can be written in terms of the following integrals:^ 



o-«i/4T 

20^ 



-u\/4t 



+ 



1-4^/ 



u+ + 4r 



(7.1) 



and 



Q{u,u',t) 



dv ipl{u)ipy{u') 







-ui/4T 



/ 



4r 



+ f 

+ 4t 



n_ + 4r 



4t 



/ 



At 



+ 4t 



(7.2) 



^Integrals P and Q may be calculated using the method of Fourier transforms combined with the convolution 
theorem. 
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2 

Here we have defined u± = u it u' and f{x) = erfc(x), with the convention erfc(x) 
1 — (2/^7r) exp{—t'^)dt. We now have 



Do 



P{u, u',n)Q{u, u",T2) + Q{u, u', n)P{u, U",T2) + 0{t^). (7.3) 



4VA 

The diffusion kernel is then 

(/ / \ \ 2— 2/-\/A 

(7.4) 

If we define 



P{u,Ti) = / du'MuTPiu,u',n), (7.5) 
Jo 

POO 

Q{u,Ti) = / du'(l)o{u'fQ{u,u,Ti), (7.6) 



then the double Pomeron source function takes the simple form 

4 /\ / 2 \ 2-2/VA 

« -jr- ("""^'1^ J "'"^"^ ■"2) + ^(^' ^2)] . (7.7) 

By evaluating P and Q we can obtain an explicit formula for TZ. Since d^Q — + P = 0, 
we need only determine Q, and use P = AQ — d^Q- For the hard- wall model the hadron 
wave functions are Bessel functions, and we can obtain an exact expression for Q, given in 
appendix 0. 

It is useful to rewrite the expressions above in terms not of n and T2 but ti±T2. Eq. ( ^.20 ) 
implies that ti + r2 = r + r_L ~ r, since t± is very small. In the double diffractive limit, the 
difference T2 — ti is proportional to the rapidity y of the glueball, defined in ( 3.11 ), 

y = VX{T2-n), (7.8) 

where we used \ti\, \t2\ <C si, S2- Since Tj is positive the range of allowed rapidity is finite but 
very large, — r\/A < y < t\/A. 

By construction TZ is symmetric in ti and T2, so it must be an even function of rapidity. 
Meanwhile the rapidity y can only appear in the ratio y/\/A. This means that the shape of 
'R{u) is nearly constant for rapidities y <C ^/X- This follows essentially from the fact that 
the gravitational scattering amplitude to which our calculation reduces in the A ^ oo limit 
is rapidity-independent. 

Let us now discuss the main features of the double Pomeron source function. TZ{u) 
vanishes in the UV as u — > oo. For large u, the Gaussians from Q and P (see eqs. (|C.4D and 
( p.5| )) dominate giving 



TZ{u) ~ exp 



-^|- + - 

Tl T2 



-TV? 
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Many aspects of TZ{u) that are true throughout the Regge regime can be understood 
qualitatively by considering the limit of large diffusion times t,ti,T2 ^ 1, and focusing 
on the behavior at small u <ti Ti. (The behavior at larger u is typically irrelevant for our 
computation, because of the exponential fall-off of the glueball wave function at large u.) In 
this limit, we have from eq. ( p.ll ) that the u and Ti behavior of TZ{u) is given by 



2 



(1 + 2n) 
(riT2)'^/^ 



u2 / 1 1 



4 \Ti ^ T2 



(7.10) 



The r- behavior was explored in ref. [5]. At U = it arises from the reflection off the hard- 
wall barrier. For modes with slow variation in u, which dominate at large Ti, the boundary 
condition required by energy-momentum conservation forces the incoming and outgoing waves 

— 1/2 

to interfere destructively, cancelling the leading behavior expected in diffusion. This 

destructive interference persists for ti <^ — A^, as the exponential behavior of the effective 
potential acts as a Dirichlet boundary condition. (Note however that at extremely large 

Ti the inverse square-root behavior is not cancelled.) 

As a function of u, TZ will have a peak at Umax = 

i[v^l-F64rir2/r - 1]. When the 
glueball has rapidity y = 0, the peak location scales as Umax ~ V^, so that the height of the 
peak scales as 7^(umax) ~ t~'^- Meanwhile, at the hard- wall, 7^(0) ~ t~^. The peak-to- wall 
ratio is 

At non-zero rapidity, keeping the leading correction \y\ <^ t\/A, the peak moves to Umax ~ 
Vr[l - (yV2T^A)] and the ratio (|7llD scales as fr[l - {y^/2T^X)]. 



We learn two important facts. At fixed central rapidity, and as r increases, the maximum 
of TZ{u) moves away from the confinement region near u = 0. This means that its computation 
becomes increasingly insensitive to the details of confinement. Furthermore, the ratio of the 
maximum of TZ{u) to its value in the confining region becomes large. This means that the 
confining region plays a smaller and smaller role in the computation of the source function. 
Unfortunately, we currently see no argument that the insensitivity of this computation to 
the confining regime should apply generally to most or all large-A gauge theories. Instead, it 
appears to be a special feature of the hard-wall model. 

In Figure ^ we plot a dimensionless version of the double Pomeron source function: 

4 / 2 \ -2+2/v^ / /7j \ -4 

a^u) = ^ ( /2io^] e-4«/v/A ^ (7_;L2) 



We have divided out the leading dependence on s and the dependence on the hadron wave 
function normalizations. We have also divided out a factor of e^"/^, which might seem odd 
since this has explicit u dependence. However, this factor is 1 at small u (since A ^ 1) 
and irrelevant at large u (since the glueball wave function falls exponentially at large u). 
The definition ( 7.12| ) is convenient in that it makes the function TZ positive-definite and 



independent of A, except through the dependence on the rapidity. 
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Figure ^ shows TZ at three different values of y; the left plot is for r = 1 and the right 
plot for T = 10. The peak, whose position in u scales as ^/T, and the ensuing Gaussian fall-off 



(7.9) beyond the peak, are clearly visible. As r increases, the height of the peak decreases, 
consistent with the (tiT2)~'^/^ scaling of the amplitude. However, the relative height of the 
peak compared to the value of TZ{0) increases. Finally, as we change y at fixed r, the shape 
of Tl{u) changes little for \y\ <^ r-v/A/2. 
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Figure 2: The dimensionless double Pomeron source function TZ, defined in eq. (7.12), plotted versus 
the AdS radial coordinate u. Each plot is at a fixed value of r and exhibits TZ at different rapidities. 
The left plot has r = 1 and shows curves for ?/ = (solid), 0.5VA (dashed), and 0.9a/A (dot-dashed). 
The right plot has t = 10 and shows curves for y = (solid), 5\/A (dashed), and S^A (dot-dashed). 
For the external hadrons, we have chosen Aq — 3, ttiq/A — 2.405.... Since A ^ 1, there is a large 
range in rapidity over which TZ is nearly rapidity-independent. 



7.2 Finite ti,t2 and large r^: absence of a Regge peak and povi^er-lavif behavior 

The scattering amplitude for our diffractive process is peaked at tj = (within the physical 
regime). In classic Regge physics, one might expect an exponential fall-off with negative tj. 
We will see in the next section that this is true, but only at small r^. In this section we will 
show that at large Tj, and for negative tj, the amplitude decreases like a power of tj, as in [2]. 

First, we show Regge behavior is absent at small negative tj. Starting with ( ^.121) , we 
calculate the amplitude as in section 7.1. In particular, we must compute integrals similar 



to those in (7T) and (7^), where there is now additional dependence on the U. We can 



only compute P and Q analytically when r, ti , T2 are large enough that the Gaussian factors 



rapidly cut off the integrals over z^. The eigenfunctions of (5.9) for all values of ti in the 



hard-wall model are given in appendix ^ We expand them to lowest order in v, to find 

HipMpie-"')] X 

[Koip^e~''') + Hipi)Ioip^e-^')] + 0(z/ 



iPUh,u)MU,u) = -u^[Ko{pie-'') + HipMpie-"')] x (7.13) 

TT 



where pi = y/—ti/A and 



-2Koipi) + piKM 
2/o(Pj) + PihiPi) 

H{pi) vanishes exponentially as pi approaches infinity; it diverges logarithmically when pi 
goes to 0, though V'!/(0, u) is finite. 
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We approximate (|5.19| ) as 

Doiu,u',u") « 2VX j di^Mh,u)i^t{ti,u') j duj^^{t2,uW^{t2,u")e-^^''^-^^'^^ (7.15) 

where we have taken t±_ to be neghgible and kept the v and w dependence only in the 
exponentials where they are multiplied by the large numbers Tj. The kernel is now easily 
computed, 

Do{u,u',u") « ^ [i^o(pie-") + g(pi)/o(pie-")][Ko(pie-"^) + H {pMPi^"'' )] 
x[Ko(/>2e-") +//(p2)/o(/>2e-")][i^o(/>2e-"") + F(p2)/o(p2e-"")]. (7.16) 

Notice that the dependence on u, u' , u" completely factorizes. Also, the dependence on r 
and Tj completely factorizes from the dependence on tj. This second factorization implies 
that there is no Regge behavior s"^*'"* ~ e"*-**^"^* in the large r and Tj limit. Looking back, we 
see this factorization of the ti and Tj dependence stems from the factorization of the v and pi 



dependence in eq. ( 7.13 ). As long as ^*(ti,u)ipu{ti, u') has a power series expansion in v near 
= 0, this factorization, and the corresponding loss of Regge behavior, will always occur at 
large Tj, for any large-A theory. 

Before continuing we should note a subtlety. The Dq kernel appearing in the above 
expression is unbounded as any of n, n', u" increase, which might seem problematic. However, 
in deriving this expression we used an approximation which breaks down at large n, u' or 
u" . One can see explicitly that when ti = 0, the corresponding expression ( [7. 3D is strongly 
damped at large u, u' , u" by Gaussian factors. Here, our use of an expansion in v, assuming 
the Ti are large, is essentially (after integrating over u) an expansion in v? jr. This is not 
valid where the Gaussian factors are important, so in this regime has no large- ti cutoff. 
However, the hadron and glueball wave functions, which are integrated against to obtain 
the source function and the S-matrix, have their own exponential tails in u. These tails cut 
off the integrals at ti ~ 1, long before our approximation breaks down. We therefore proceed 
without concern. 

Next, we show that the scattering amplitude falls off as a power law in tj/A^, another 
universal phenomenon. The presence of falling powers is easy to understand. At large negative 
values of tj, the effective potential in the Hamiltonian (|5.9| ) develops an exponential barrier 
that forces the corresponding analogue particle away from the wall. In the calculation of the 
scattering amplitude, the eigenfunctions of this Hamiltonian are integrated against the wave 
functions of the glueball and external hadrons, which have falling power-law tails in z at small 
z (exponentially falling tails in u at large u). For large values of tj, the eigenfunctions lack 
support near the wall, and the computation is dominated by the power-law tails of the wave 
functions. Once this is true, the entire computation scales with tj, giving the amplitude a 
power-law dependence on ti. 

To demonstrate this explicitly, we insert ( 7.16| ) into ( |5.16| ). We are left with two factorized 



integrals over n' and ti", which we will denote /' and /". In the limit of large /Oj, H{pi) vanishes. 
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so the integral I' becomes 

/ du'M^'fKo{pie''^)= —Mz'?Ko{^/\hW)■ (7-17) 
Jo Jo ^ 



The integral /" is identical with t2 replacing ti. Since the function Kq{x) is exponentially 
damped for large x, the integrand only has support for z' <^ l/y^\ti\. In this region, the 
external hadron wave function (f)o{z') has a power-law tail: 

00(^0^ /vf'la/^ ivA iz'«m,'^zo). (7.18) 

VVolwi?^/^ r(Ao - 1) Jao-2("^o/A) 

(Recall that, for low-lying external hadron states, niQ ~ A = 1/zq.) With this approximation, 
/' can be evaluated by extending the range of integration over z' from zero to infinity. Inserting 
the integrals I' and I" into the definition of 'JZ{u), eq. ( 5.16| ), gives 



x(Ao,mo) = ( . (7.20) 



where 

x(Ao,mo) = ( /A\ 2 

This expression is valid in the region of our approximation, u <^ y^, as noted above. Note the 
almost-quadratic growth of the amplitude with s, the power-law dependence on the ti, and the 
iTiT2)~^^'^ factor, which provides subleading s dependence and limited rapidity dependence. 
The two Bessel functions determine the shape of the function in u and contribute a subleading 
dependence on the tj. Since pi = -y/jtij/A ^ 1, these functions assure an exponential cutoff of 
the source function near the wall at u = 0. The growth of the source function as n ^ oo is cut 
off by the breakdown in our approximation. But before this happens, the growth is more than 
compensated by a falling glueball wave function {(j)5{z) ~ z^) in our current computation. 

We can now compute glueball production at large ti, using the above Pomeron source 
function. The amplitude involves the integral of the glueball wave function against the 
Pomeron source function. The factors of KQ{y^\ti\z) are damped exponentially near the 
wall, so the glueball wave function can be approximated by its power-law tail. Also, the 
factor e^"/^ can be replaced by 1. 

The scattering amplitude (5J) obtained from (|5.17] ) is 



where 



2x(Ao,mo) /m\2 /m^X^" A^ 



J2{m/K) \AJ \t1t2J (ti-i2) 



-l + l^ln(^ 

2ti-t2 \t2 



(7.22) 
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The amplitude is well-behaved as t2 ^ ti, since 

x(Ao,mo) fm^^ m^^"A^ 

For equal momentum transfers, the scaling of our amplitude in terms of gauge theory variables 
is simply 

Although we are not quite in the same limit, we may compare ( |7.24| ) with the high energy, 
fixed-angle scattering result of ref. [2]. At large s and ti but fixed s/ti, eq. (14) of ref. [2] 
indicates that in 2 — > 3 scattering of four objects created by operators of dimension Aq and 
one of dimension A = 4, the amplitude would scale as 

where p is a characteristic momentum scale for the 2 — > 3 scattering process. General require- 
ments assure the overall power of gs ~ is the same in both cases as well as the 1/A sitting 
in front. For low-lying hadrons m ~ rriQ ~ A. Our amplitude scales as exhibiting 



a power law suppression in momentum with the same overall exponent as ( 7.25| ); both scale 



as (A/p)^^o. However, the dependence on \/A = B? /a' in the two expressions is different. 
This is because the dynamics of string theory itself — the exponential suppression of hard 
scattering — cuts off the amplitude in [2] , which introduces factors of a' . Here the cutoff on 
the amplitude occurs through the Kq Bessel functions above, through the large momentum 
transfers ti and t2- These are independent of a' . 

7.3 Small ti, t2 and small t^: classic Regge phenomenology 

In the classic Regge regime of s ^ |t| ~ A^, experiments show hadronic amplitudes exhibit 
Regge behavior, S ~ s"^*) with a{t) approximately linear in t. We have just seen that for 
large r, this behavior is absent. What has happened? In [2] it was argued that for fixed t the 
classic Regge behavior of the amplitude is a transient phenomenon present only for a certain 
range of s. When s is sufficiently large, a transition takes place and power-law behavior is 
restored. Here we will see a similar phenomenon. In this section, we show that for small 
enough ti and Tj, one finds approximate Regge behavior, as a transient effect that gives way, 
for larger Ti, to the power-law behavior seen in the previous section. More specifically, we 
will work in the regime 

^ < Ti < 1 and ^ ~ 1 . (7.26) 
The lower bound on is the double Regge limit we have taken from the beginning. The 



upper bound on Tj keeps the diffusion times short enough that the potential ( 5.10 ) can be 



- 20 - 



approximated by its first few terms in a Taylor series expansion. The second relation keeps the 
product Ti\ti\ large enough that we can distinguish between exponential and linear behavior 
of the amplitude as a function of ti. We will argue that the presence of this window, where 
Regge scaling can be seen, is universal and is not sensitive to the details of the hard-wall 
model. 

The generality of the phenomenon can be seen from the following argument. The source 
function TZ{u) involves the combination of two Pomerons, one from each of the external 
hadrons, 

7^(n) ~ {u\e-^^^^\<t>l){u\e-"^^^\c^l) , (7.27) 
each of which requires a diffusion-kernel calculation of the matrix element 

{u\e-"^^^(t)l) = [ du'{u\e-"'^^u'){u'\cbl) (7.28) 



where Hi is the Hamiltonian ( |5.9| ) and |0q) is the state whose wave function is (poiu) in the 
position basis. That TZ{u) is given by (|7.27| ) can be seen from examination of ( 5.16 ) and 



For small Tj and ti we may evaluate the matrix element ( [7.28|) semiclassically, writing it 
as a path integral over all paths between u and u' , and expanding around the classical path 
of minimal Euclidean action. Defining the average position u = {u + u')/2 we may rewrite 
the potential as 

Vi{u') ~ 4 - zgi,e-2"[l - 2{u' -u) + 2{u' - uf + ...] . (7.29) 
Working to linear order in the potential, the matrix element is then 

(7x|e-^»"»|0g) = j du'exp [zitine-^^] e-^^^c^liu') J Vy e- fo' ds[ly^+lp{y-u)] _ ^^^gg) 

Here /? = 4zQtje~^", and the boundary conditions on the path integral are that y{0) = u' and 
y(rj) = u. The classical solution is yci{s) = ^Ps'^ + vos + u' where vq = —\(5Ti + '^^-^. Letting 
y = Vd + X, the path integral becomes^ 



du'exp [z^tine-^''] e"^"'+^^ ~<Pliu') / Vx e- ^o' i^'^s^ (73^) 

where the boundary conditions on x are x(0) = x(rj) = 0. The remaining path integral 
is proportional to Xj ^/Ti but is independent of ti and plays no further role. The Gaussian 
factor in the integral, along with the constraint on Tj ( 7.26| ), constrains v! to be near u. 



^Keeping the quadratic term in the potential Yi(v!') wiU produce contributions to the transition am- 
plitude that are subleading in Ti. One might worry that it is inconsistent to consider the quadratic i? 
fluctuations without the fluctuations. Keeping the fluctuations, the quantum determinant scales as 
|/3i |'"^*csch(2\/2|/3i Iri)"*"^^ which is indeed independent of (ii for small Ti (see for example [16]). 
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Expanding the rest of the integrand around u' = u, performing the Gaussian integral, and 
dropping terms suppressed exphcitly by powers of Tj, we obtain, after substituting the exphcit 
form for (j)o{u), 

{u\e-^^^^\^l) ~ e-4"-^»'=-'"jAo-2([mo/A]e--)2erfc +0(^), (7.32) 



where bi = —ZQtiTi > 0. The Pomeron source function TZ{u) is then proportional to (7.32) for 
i = 1 times a similar expression for i = 2. 

Finally, to obtain the glueball production amplitude, there is still the integral over u to 



perform, (5.17), of the Pomeron source function against a glueball wave function. That this 
integral leads to approximate Regge behavior, with a Regge slope that slowly varies with tj, 
can be seen as follows. The complementary error functions can be replaced by 2 everywhere 
in the range of integration except for u < ^/Ti, a region which gives a subleading effect in Tj. 
The integral ( 5.17] ) is dominated by the exponential exp(— 6je~^") and the power law tails that 



come from expanding the Bessel functions at large u. Thus a reasonably good approximation 
to the ti dependence of the amplitude is 



oo roo 1 







duMu)n{u)^ I due-"^ = -6-*^ (r(M) -r(M,6)) , (7.33) 



where r(a, z) = t"" ^dt is the incomplete gamma function, b = 61+62, and M = 2Ao+2. 
Expanding this result at large M leads to approximate Regge behavior, e^^^^^^, where we 



recall that bi oc — ZQtjTj. 



In sum, we find approximately exponential behavior in the window ( 7.26| ), and therefore 



Regge behavior of the amplitude, with a Regge slope of order 1/\/AA^. Once Tj or ti becomes 
large, however, this argument breaks down, Regge behavior is lost, and the physics enters the 



regimes discussed in sections 7.1 and 7.2 



8. Discussion 



Using the hard- wall model for gauge/string duality, we studied glueball production from 
hadron scattering in the double diffractive limit. String theory provides a new window into 
this regime, which is not accessible using perturbative QCD. A principal objective of this paper 
was to set up a formalism that can be applied to a number of 2 ^ 3 processes that are more 
experimentally relevant than glueball production. These could include heavy quarkonium 
production and Higgs production. 

Our study is therefore largely of technical interest, though we also uncovered some notable 
physical phenomena. We wrote our amplitudes in terms of a source function TZ{u), which 
was independent of the glueball state and characterized the two Pomerons which fuse to 
make the glueball. The amplitude is given by integrating this function against the glueball 
wave function. Applying this formalism, we saw that Regge phenomena are only present at 
relatively small s, si,S2, and are lost as these quantities become large. We found that at large 
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s, Si and small ti, the scattering amplitude becomes rapidity- independent, a fact which arises 
from the rapidity-independence of the corresponding tree-level dual gravitational amplitude. 
We also found the expected power-law fall-off of the amplitudes at large momentum transfers 
ti, where scaling behavior is expected, as in [2]. 

A natural next step would be to apply this formalism to heavy quarkonium production. 
Quarkonium in AdS / CFT can be modeled by adding a D-brane to the AdS cavity that fills 
AdS from the boundary to some minimal radius ro = / zq [17]. The mass of the heavy quark 
scales with ro, and the quarkonium supergravity wave function must live on the D-brane. In 
principle one ought to consider the double Regge limit of the five-point disk amplitude for 
string theory in flat space, with four closed string insertions, and one open string insertion 
representing the heavy meson. However, most of the important physical effects may be 
captured simply by using the fact that the quarkonium wave function has support only in 
the region r > tq {z < zq) ignoring any further details. Even with such a simple model, it 
might be possible to make a prediction for the relative rates of bottomonium and charmonium 
production in the double diffractive limit. 

Double diffractive Higgs production might be a clean way to observe the Higgs boson 
(or other scalars) at the LHC. There are competing field theory models, some inspired by 
flat-space string theory, which predict potentially observable cross-sections for exclusive Higgs 
production at the LHC. See ref. [8] for a brief review. It is possible that gauge/string duality 
might clarify some of the approximations made in these models. More concretely, we would 
proceed as follows. The Higgs, which couples to F^yF^^^ in the standard model through a 
dimension-five operator, should be treated similarly to the scalar glueball in our discussion 
above, except for one crucial difference: we should replace the normalizable dilaton mode in 
supergravity, representing the glueball, with its non-normalizable counterpart, representing 
the Higgs boson. This non-normalizable mode, at timelike momentum p^^ with —p^ = rnjj, is 
highly oscillatory. The technical challenge that lies ahead is then to compute the fusion of two 
Pomerons into such a state. We expect that this challenge can be met, and that gauge/string 
duality will soon contribute to the debate over diffractive Higgs production. 

Acknowledgments 

We would like to thank Andreas Karch, Richard Brower and Chung-I Tan for discussions. 

The work of M. J.S. was supported in part by the U.S. Department of Energy under Grant No. 
DE-FG02-96ER40956. The work of C.P.H. was supported in part by the National Science 
Foundation under Grant No. PHY-0243680. The work of S.P. and E.G.T. was supported in 
part by a Royalty Research Fund award from the University of Washington. 

A. Bosonic tachyon/dilaton amplitude 

In this appendix, we calculate the double diffractive limit of the four tachyon, one dilaton, tree 
level, closed string amplitude in bosonic string theory in twenty-six dimensional Minkowski 
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spacetime. We let the dilaton correspond to the fifth particle in Figure |I|. Following closely 
an approach described in [12], we are able to express the result both as a power series in a'm^ 
and an asymptotic series in l/a'm]_. The power series result is important for our discussion 
of glueball production in the paper; for glueball production, a'm^ is effectively very small. 

We use the conventions in Polchinski's textbooks [18]. In section |5[ string theory momenta 
have a~ to distinguish them from four-dimensional gauge theory momenta. We will not need 
to make that distinction here as we do not work with the four-dimensional gauge theory 
momenta at all. Unlike the conventions in the main body of the paper, we take all of the 
momenta ki to be ingoing. 

Take four tachyon vertex operators Tj = gc :e^^^'^:, i = 1, ... ,4 with kf = 4:/a', and a 
dilaton vertex operator D = g'^-f^v idX'^dX'^e^^^'^: with /c| = 0. The dilaton is a massless 
particle that travels on lightlike geodesies so the symmetric tensor /^j^ must be transverse 
to its momentum, k'^f^y = 0. Given another lightlike vector k such that /cs • A; 7^ 0, we can 
satisfy the constraint by taking f^y = rj^y — {k^^j^ky + kfj_k^y)/k^ ■ k. The correlation function 
of these vertex operators on the Riemann sphere is 



I 'V . . I ^ ki'kj 



\i=l 




We have used the shorthand notation Zij = Zi 



The tree-level S-matrix for the scattering of four tachyons and a dilaton is obtained 
by integrating the correlator ( |A.1| ) over all possible worldsheet coordinates for the operator 
insertions, weighted by a topological factor. 



-2m/" (fzi Agi^ost (T[Ti(z^, Zi)D{zr„Z5)) . (A.2) 

:!li^cu{oo} \._, / 



The path integral over ghost fields contributes a Jacobian Aghost = C^i^'^i^a — ^2)5^(2;, — 
— z'^)\zabZbcZca\'^ that fixes the residual PSL(2, C) symmetry left over from conformal 
gauge-fixing. This gives us freedom to fix three vertex operators to arbitrary positions. A 
convenient choice is zi = 0, Z4 = 1,2:5 = 00, because that makes 1-2142:45 251!^ = |z5|^ and 
ni<5 kisi" ^''^^ — 1- Since ( |A.1| ) comes with a factor of |25|~^, the S-matrix is finite.^ 

Tse = C / d^Z2(fz^ J] \zijr' Uu{k2Z2 + k^z^ + kiY{k2Z2 + A:3Z3 + k^Y . (A.3) 
i<«<i<4 

We have from ref. [18] that e~^'^C_^C^2 = '^t^ jo! g^ and = 2gc/a' , which leads to the result 
that C = —Airg^. In these conventions, the coupling 2TTgc = k, where k is the gravitational 



®Note, ii z — X + iy, we define the measure factor d?z = 2dx dy. 
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coupling constant that appears in front of the Einstein-Hilbert action as {2k?)^^. Thus, we 
find that gc ~ a'^'^~'^^^'^gs where gg is the dimensionless string coupling (which is related to 
the vev of the dilaton) and d is the number of dimensions in which we work. 

Make the conformal transformation u = z^jzj, and t; = (23 — l)/(z2 — 1) to get 

T26 = C I d2nd27;|nr"'*^/2-4|^|-a't2/2-4|-^_^|-a'si/2-4|^_^|-a's2/2-4|-^_^^|a'(-s+si+S2)^^ 

[k2uil -v) + ksil -v) + kiil - uv)f [c.c.]\ (A.4) 



Now take the double Regge limit given in ( |3.6D . The dominant integration regions are near 
the origins of the u and v planes. In order to demonstrate this we should first discuss the 
issue of convergence. For fixed v and for real Mandelstam variables with physical signs, there 
are clearly three special points in the finite n-plane where the integrand becomes singular, 
u = 0,l,v^^. Of these, the origin is benign as long as we take a'ti < —4 so that the 
singularity is integrable. The other singularities at u = l,v~^ are not integrable for positive 
si and s. One way to avoid this difficulty is to choose pure imaginary values for the s- 
type Mandelstam variables and analytically continue the result to the physical domain. In 
this scheme |1 — n|~"'''i/^, |1 — 't;|~"'*2 ^nd |1 — ot|"'(~*+**i+'*2)/2 are just phases. Now if 
a'ti > —8, then the integrand vanishes sufficiently fast at infinity for the u-plane integral to 
converge. Thus, the entire integral may be defined by continuation from —4 > a'ti > —8 
and pure imaginary s,si,S2- We observe that this range of momentum transfers implies that 
|^|-Q!'ti/2-4 |^|-Q't2/2-4 always singular near the origins of the u and v planes. 

For fixed v, the integral in u is dominated by a saddle point at u ~ 0(l/si), which 
approaches the origin in the Regge limit. Therefore, |1 — n|~"'*^/^~^ ^a'si{u+u)/4 ^ Similarly, 
for fixed u, the integral in v is dominated by a saddle point at v ~ 0(l/s2). It follows that 
|1 _ ^;|-"'^2/2-4 ^ ^a's2{v+v)/4 and 1 1 - ^it,|"'(-^+^i+«2)/2 ^ ^a' s(uv+uv)/4 _ Therefore, 

T2&-C j d\d2r;|u|-"'*i/2-4|^|-a't2/2-4g^p \^a' siRe{u) /2 + a' S2Re{v) /2 + a' sRfi{uv) /2\ 
xUu [k2u{l -v) + ksil -v) + k^{l - uv)f [c.c.]" . (A.5) 
We still need to understand the consequences that the Regge limit has on the second line 



of (A.5). It suffices to consider only the r/^,/ part of the tensor since the portion containing 
k corresponds to a longitudinal polarization that must decouple in physical processes [19]. 
Moreover, k is arbitrary and the final amplitude cannot depend on it. So 

Vfiu [k2u{l -v)+ ksil -v)+ A;4(l - uv)]^ [c.c.]" 
= {k3 + k^f + k2 ■ {ks + k4){u + u) - k?, ■ (fcg + ki){v + v) 

+ (terms with two or more factors of u,u,v,v). (^-6) 



Scaling u — > u/si and v — > v/s2 in ( A.5 ), it is clear that only the first three terms on the 



right hand side of (|A.6D have the possibility of not being suppressed by a large energy.'' 



careful treatment of this scaling is given later in this appendix. 
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Since ■ (k^ + k^) = —^2/2, the third term is actuahy suppressed. Using k2 • (k^, + k^) = 
(^2 — si — 4/a')/2, we may replace ( |A.6| ) by —t2 — si{u + u)/2. 
For convenience we define a prototypical integral 

/(a, a, b,b) = j (fu (fv u^lfv^lfi exp [a s{R,e{u) /2 + a's2Re(f ) /2 + as Re{uv) /2] . (A.7) 
With a = —a'ti/A — 2 and b = — a't2/4 — 2, the amplitude reads 

T2&/C « -t2l{a, a, fc, - y [/(a + 1, a, ^ + /(a, a + 1, 6, (A.8) 

726 is a function of five variables s, si, 52, ti, ^2- It is still not well-defined for physical scattering 
(real s > si + S2 > and real ^1,^2 < 0) since I diverges. For a proper definition, we extend 
s,si,S2 to the complex numbers and follow a technique due to Lipatov [12]: decompose I 
over different regions of the u and v planes, and for each region choose the phases of s, si, S2 
such that they lie on the real axes with signs that ensure the convergence of the integral. 726 
is then defined by analytic continuation of s, si, S2 to the positive real axes. 

We now carry out Lipatov's procedure for / as given in ( |A.7] ). It is convenient to let 
a' = 4. At any fixed value of {v,v) the integrand has a saddle point at {u,u) = (— a/(si + 
sv),—a/{si + sv)). Similarly, for any fixed value of the integrand has a saddle point 

at {v,v) = (— 6/(s2 + su),—b/{s2 + su)). These saddle points lie close to the origin of the 
complex plane. 

We begin by evaluating ( [A.7] ) directly which yields an asymptotic series in (a'm^)~^. 
Such an expansion is useful in the regime a'm^ ^ 1, but the glueball production process 
investigated in the paper involves the opposite regime a'rn^ <^ 1. Thus, we follow by changing 
the variables of integration in ( [A.7| ) to yield an integral over uv and a power series in a'm^. 

A.l Large a'm'] 



We divide up ( A.7 ) into four pieces, depending on whether u and v have positive or negative 
real parts, 

/ = ^/'^".'^" , (A.9) 

where we have introduced the notation ax = sgnRe(x). In each domain, we choose si and 
S2 to have opposite signs from Re(M) and Re(u), respectively. In this way, the exponential in 
( [A.7] ) damps, ensuring that the integral converges. 

For /'^'"'^"j we choose si = |s]^|e*'^(°""+^)/^ and S2 = \s2\e^'^^°'"'^^^/'^ and change variables to 
w = e~^^siu, W = e~^"'^^ siu, z = e~^'^S2V, and z = e~^'^^'^ S2V. These peculiar transformations 
guarantee that w is the complex conjugate of w and that 'z is the complex conjugate of z. 
Then we find 

I<ru,o, = (e-*'^si)-"-i(e-*'""^si)-^-^(e-^'^S2)"^"^(e"^'""^S2)"^"Vi (A.IO) 

where 



Jl= j d^wd^zw''ufz''z^exp 

Re{w)>0 
Re{z)>0 



-{w + w) — (z + z) -\ {wz + Wz) 

S1S2 



(A.ll) 
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where® 



Now Jl depends only on the combination s/siS2- Restoring a', this is the dimensionless 
ratio 4/a'm^. In the hmit where a'm\ ^ 1, we may expand the exponentials in ( |A.ll ) as a 
double power series, 

{a'm? /4)~"-~™ _ 
Jl= y ^i-n K(a + n,a + m)K{b + n,b + m) (A.12) 

n,m>0 

Kix,y)^ [ d^u,w-lBye-^''<^^ = 2n^^^^^^%P^. (A.13) 

jRe(w)>Q ^{-X)T{-y) 

The last equality holds only for —2 < Re(x + y) < —1, but the result can be analytically 
continued outside of this range. This is necessary in order to evaluate K for physical values 
of a'ti as well as nonzero n and m. 

A.2 Small a'm^ 

As we are interested in in the opposite regime a'mP: <C 1, we shall not process the large 



a'm^ result further. Instead, we make a change of variables in ( A. 7 ) in order to be able 
to derive a power series in a'm^. We begin by switching the integral over u in (A. 7) to an 
integral over uv. We will see that this change of variables misses the saddle point at small 
\v\. Thus, we will eventually need to add a second contribution where we switch instead the 
integral over v to an integral over uv. 



v-plane 




u-plane 



Figure 3: Schematic representation of the dominant integration domains for /. The vertical strip 
intersects the u-plane in a small disc surrounding the origin (m, u) = (0,0), while the horizontal strip 
intersects the i;-plane in a small disc surrounding the origin {v,v) — (0, 0). 



We will call Ivcrt the integral (A. 7) in which we have replaced an integral over u by an 



integral over uv. The subscript "vert" indicates that we are picking up the saddle point in 



*After changing to polar coordinates and rescaling, the integral K factors into the product of two one- 
dimensional integrals and may be readily evaluated. 
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the vertical strip in Figure ^ In switching the integral over u to an integral over uv, we can 
divide up the domain of integration into the four regions where Re(?;) and Re(nu) are either 
positive or negative, 

/.ert = ^C;r . (A.14) 

We begin by considering the domain where Re(v) and Re(tif ) are both positive. This neces- 
sitates choosing S2 = \s2\e^'^ and s = |s|e"^ so that converges. Making the change of 
variables w = e~^'^suv and w = e~^'^suv, we find 



vert 



(A.15) 



in „\ — 2—a~a 



(fwcfv w^llfv''-''-^lf'-^-^ exp 



Sl / W W , , , _x 

h — + S2{V ^V)-(W^W) 

S \ V V ' 



Re(w))>0 
Re(D)>0 



Changing variables to z = e ^^S2V and z = e ™ S2V gives 
where 

'"siS2 I w w 



Js 



(fw(fz w^Ufz^-"-^^-'^-^ exp 



+ 
z z 



{z + z) — {w + w) 



(A.16) 



(A.17) 



Re(«))>0 
Re(z)>0 



Next consider the integral Z^^- Now we must choose S2 = |s2|e*'^ and s = |s| so that the 
integral can converge. Change variables to w = suv and w = suv, then to z = e~^'^S2V and 
z = e~^'^S2V, and finally rotate w — > e~'^'^w and w — > e^'^w. This yields 



^vcrt ~ 



— 2—a—a/—i-K \a+a—b—biiT(a—a) 



{e -82)"""' " "e"-'" "'Js- 
Thus, for the two terms in ( [A.14| ) with Re(t') > we have found 



vert * vert 



iyr j,^ — 2— a— a _j_ ^iTT{a—a)^—2—a—a 



{e''^S2T^~^-^-^Js- 



a+a—b~b 



Js- 



A similar analysis for the four terms with Re(f ) < gives 

J — h ^ _ ^iiT[{b-b)-{a-a)] ^^-iiT ^^~2~a-a _|_ ^iiT{b-b) ^~2~a~a 

Combining (|A.19| ) and ( |A.20D , 

f —in \ — 2—a~d _|_ i7v{a—a) —2—a—d ftD~'^'^r. ^a^d—b—b _[_ i7v[(b—b) — {a—d)] a-\-d—b—b 



(A.18) 
(A.19) 
(A.20) 



'vert 



Js- 
(A.21) 

Like Jl, Js depends only on the combination S1S2/S. Restoring a', this is the dimen- 
sionless ratio a'm^ /4. As such, we may expand the exponentials in ( [A.17| ) as a double power 
series, 

Js= ^1;^ K{a + n,a + m)K{b-a-l-n,b-a-l-m) . (A.22) 



n,m>0 



n!m! 
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Simplifying and using the identity 7rcsc(7rz) = T(z)T{l — z) we obtain 



Js 



47r^ sin(7ra) sin [7r(6 — a)] 
sin [7r(a + a)] sin [7r(a + a — 6 — 5)] 



E 

n,m>0 



(a'm2_/4)"-' 



r(i + a + n)r(6- 



n 



n\m\ 



r(— a — m)r(— 6 + a + 1 + m) 

(A.23) 

The form of ( |A.17 ) along with our power series method of evaluation makes clear why 
/vert is not equal to /. We have assumed not just that a'm\ <C 1 but that a'm\_w/ z ^ 1 in 
the domain contributing to the integral. However, this condition cuts out the small z region, 
which in turn corresponds to the small v region, which is where the other saddle point of / lay. 



To cure this problem, we will instead switch in (A. 7) the integral over v for an integral over 
uv. This second term /horz will contain the saddle point at small v but will miss the small u 
saddle point. The subscript "horz" indicates that we are now picking up contributions from 
the horizontal strip in Figure |3|. The complete result for / is given by the sum /vert + /horz- 
Repeating steps (|A.14| ) to ( |A.23|) for 

-^horz is straightforward. In fact, /horz can be obtained 
from /vert simply by exchanging a ^ b, a <^ b, S2 <-> si. This completes the solution of 726. 
Suppose 



a'm^i < 1 



(A.24) 



so that it suffices to keep only the n = m = term of ( |A.23| ). Then /horz(a + 1, a, b, b) = 
/horz(a,a + 1,6, 6) = (6 - a)/horz(a, «, ^, Also, /vert (a + 1, a, 6, 6) and Inertia, a + l,b,b) 

are each proportional to (s2/s)/vert(ffli b, b). When they are combined with the factor of si 
in ( [A.8| ), they contribute at order a'm\ which we are dropping. Thus, 



/26/C~ -t2/vcrt(a,a, 6, 6) - tihorzia,a,b,b). 



(A.25) 



This is our main result. Using the identity (e ^'^z)^ + z^ = 2{e ™/'^z^^ cos(7rx/2) and restor- 
ing factors of a' gives 



To 



26 



a' 



+ e 



-i-K/2 



-i7r/2 


a's\ 




~) 




a'ti/2 


a' s\ 


-) 






a't2/2 


a' s\ 







t/2Q ^2 
4 



4 



(A.26) 



, X a'(t2-il)/2 



n(a'ti/4,a'(t2-ti)/4) 



, X a'{ti-t2)/2 



n(a't2/4,a'(ti -t2)/4) 



where 11 is defined in eq. (5.4). 

We can check that this amplitude is consistent with expectations from supergravity by 
expanding the IT functions near the graviton poles at ti = and t2 = using eq. We 
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find 



26 



a' 



+ 



+ 



4 



1 



2+a'(ti+t2)/4 



a'(ti-t2)/4 2sinh 



a'ti/4 
1 



i7r/2 








j7r/2 









2+«'t2/2 / / ^ a'(ti-t2)/2- 



(A.27) 



If a'ti and o/t2 are slightly less than 0, then the first term remains finite, whereas the second 
and third terms diverge and the corresponding Regge exponents for s will be slightly less than 
2. 



B. Eigenfunctions in the hard- wall model 

Consider the eigenvalue problem Hi'il}{u) = E'(p{u), where Hi = —d^ + 4 — ZqUc"'^'" and 
E = 4 + u'^ for u > 0. The form of the operator Hi was derived in [5]. 

As a shorthand, denote the dimcnsionless momentum transfer p = Zo\/|^i|- Then in 
terms of the variable ^ = the Schrodinger equation is the modified Bessel differential 

equation: 

eV'+eV''-(^' + (^^^)')V' = o. (B.i) 

The general solution is in terms of modified Bessel functions of the first kind, 

M^) = c {liAO + Ri^, p)i-iAO) ■ (B.2) 

The relative coefficient R{i',p) is determined by the boundary condition at the wall [5], 

m'^)k=p = ^- (B-3) 

This boundary condition follows from energy-momentum conservation. More precisely, note 
that the metric fluctuation = is pure gauge because it corresponds to a linear 
reparametrization of the background metric. In order to preserve diffeomorphism invariance in 
the bulk and correspondingly energy-momentum conservation in the boundary gauge theory, 
we must require that this pure gauge metric fluctuation satisfy the boundary conditions at 
the hard wall. 

The boundary condition yields 



_ ^hujp) + pliv-l{p) + pliu+l{p) 

^=p + (p) + pl-iu+i (p) ' 



(B.4) 
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Note that at small u, R{v, p) —1. This factor of minus one leads to destructive interference 
in ipuiu) near n = 0. In section ^, our integrals over ipy{u)dv are dominated by small v, 
and thus this destructive interference has important effects for the qualitative features of the 
scattering amplitude. 

The overall coefficient c is fixed by requiring that the eigenfunctions are delta-function- 
normalized in the coordinate n, 

oo 

dv il)l{u)ipy{u') = 5{u - u') . (B.5) 





This gives |cp = i//(2 sinh(7ri/)).^ This normalization does not fix the phase of c. We can 
choose this phase so that ipuiu) has a well-defined limit as p ^ 0: 

f p\ -i-v vViiv) 



Given (|B.4D and ([B.q ), we have solved the eigenvalue problem Hi'il){u) = E'4j{u) 



In section [7.1| , we need these eigenfunctions in the = limit. In this limit of vanishing 



momentum transfer, we Taylor expand around /3 = to get 



Mu) = f " + ^e-" ) + 0{p'). (B.7) 
V27r V u + 2i 



C. Calculation of Q 

The integral 



Q{u,Ti)= du'MurQ{u,u',n) (c.i) 







may be done by expanding the external hadron wave functions (j)o as a power series in e " 
near the UV boundary of our cavity: 



oo 



n=0 

The leading term in the expansion will capture the UV behavior of the wave function and 
higher order terms will capture the IR behavior — the details of which depend on the precise 
physics of confinement. For our hard-wall model, (j)o{u') is a Bessel function where the odd 
coefficients vanish, C2n+i = 0, while the even coefficients are given by 

mo/A * " 



, ^^^^2_i^j (-1)- ,^3. 

AVVolvKi?3/2 JA„-2(mo/A) n!r(Ao-l + n)' ^ 



®This normalization constant can be obtained by working at small momentum transfer p where the eigen- 
functions look like plane waves. 
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Inserting the power series for (t>o{u') in eq. ( p.l|) , the resulting expression for Q can be ex- 
pressed compactly as 

/"OO 

Q{u,Ti)= CnCm du' e-^' Q{u,u' , T,) . (C.4) 

n,m>0 

The integral over Q(u,u' ,Ti) was then done exactly using Mathematica 5.2: 



noo 



e 



Jo 16 X - 1 



u + 4Ti\ + 3 „ / li + 4r, 



4rj / X - 1 V ^/4Ti J X + 1 V ^/Wi 

where /'(y) = --^ + 27//(y). 

We now consider the limit of Tj 1 with n <C Tj in which case we can obtain beautiful 
formulae for Q, P, and TZ{u). In this limit, we apply the asymptotic expansion f{y) ~ 

(y^ ~ 2^" ~'~ if^ + • • to eq. ( p.5| ). We find that the integral, and hence the summand in 
eq. ( |C.4| ), scales asymptotically as r- 

rdu'e-'-'^Q{n,u',n) ^ il±^l±^ e""'/^- . (C.6) 

With these formulae and in this limit, we can understand the Tj and u dependence of Q: 

QKr.)-^i^e-"V^-C (C.7) 

where 

^_ V- 2(2 + 2Ao + n + m) 
^ (2Ao + n + m)^ 

n,m>0 

With knowledge of the c„, C can be calculated numerically and rescales Q by some overall 
constant. For example, for our hard-wall model, in the case mo/A = 2.405... and Aq = 3, we 
find that 

" A2VoVi?3 (7-270 X 10-3) (l + 2n)^-37^ " ^^'^^ 



(For the hard wall, there is in fact an analytic formula for C as a function of Ag and mg/A 
involving 2F1 functions.) 

To obtain an expression for TZ{u), recall that P = 4Q — d^Q- In the large Tj limit, 
the derivative is subleading, and we obtain P ~ 4Q. TZ{u) in the large limit can then be 
constructed from (^): 

- ^ (^^^'^^'^J e^-/v^g(n,n)Q(n,r2) . (C.IO) 
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Neglecting the e^"/^ because of the large A limit, the dependence of this expression for TZ{u) 
on u and Tj is 



n{u) 



{i + 2uy 



■ exp 



u 

'T 



1 1 

n T2 



(C.ll) 



(riT2)3/2 

In order to see how accurately approximates ( |7.7| ) , in Figure ^ we plot the dimen- 

sionless version of the double Pomeron source function, defined in ( [7.12 ), for various values 
of Ti and T2. As expected, the agreement gets better as each Tj increases. 




Figure 4: We plot TZ{u) for various values of ti and T2. The solid curve is the exact numerical value 



of TZ{u) as a function of u. The dashed curve is our approximate formula ( C.ll ). We have plotted a) 
Ti = T2 = 25, b) Ti = 10, T2 = 40, and c) ti = T2 = 5. The approximation, always too large, improves 
as Ti increases. 



We also examine the accuracy of eq. ( 7.11 ) by plotting in Figure ^ its predicted value for 
^('"max)/^(0) against the actual value obtained numerically from (|7? 




Figure 5: We plot the ratio Tl{u„^ax) /T^{0) as a function of t/2 with ti = T2. The points are the 



numerically determined values of the ratio, while the solid line is the approximation eq. (7.11). 
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